Correlations and superfluidity of a one-dimensional Bose gas in a quasiperiodic 

potential 



O 
(N 

< 
in 



I 



cd 



I 

o 
o 



cn 

in 

o 
o 



% 



Alberto Cetoli and Emil Lundh 

Department of Physics, Umed University, SE-90187 Umea, Sweden 

We consider the correlations and superfluid properties of a Bose gas in an external potential. Using 
a Bogoliubov scheme, we obtain expressions for the correlation function and the superfluid density 
in an arbitrary external potential. These expressions are applied to a one-dimensional system at 
zero temperature subject to a quasiperiodic modulation. The critical parameters for the Bose glass 
transition are obtained using two different criteria and the results are compared. The Lifshits glass 
is seen to be the limiting case for vanishing interactions. 

PACS numbers: 



I. INTRODUCTION 

Research on the problem of a Bose gas in a disor- 
dered potential gained momentuni with the realization of 
well-controlled disordered potentials for ultracold bosonic 
atoms, which inspired a surge in theoretical and experi- 
mental activity [1|-u|j ^^.d the subsequent observation of 
Anderson localization by two groups [g, 0] . The quest is 
now on to better understand the effect of disorder on an 
interacting quantum system. 

The phase diagram of an interacting Bose gas subject 
to a disordered potential was outlined in, e.g., Refs. [^ 
0, [iOj [ill • While the exact picture of the phase diagram 
depends on the employed potential, the qualitative be- 
havior appears to be the same for all kinds of disorder 
realizations. In absence of disorder, the gas is in a su- 
perfluid phase, which in a two- or three-dimensional Bose 
system at zero temperature is identical to a Bose-Einstein 
condensate (BEC), and in one dimension is a quasi- 
condensate with suppressed density fluctuations and al- 
gebraically decaying correlations. Upon raising the dis- 
ordered potential, the system can enter a new quantum 
phase: the Bose glass state, also called a fragmented 
Bose-Einstein condensate. The Bose glass state lacks 
long-range phase coherence, having exponentially decay- 
ing correlations and zero superfluid density, but it is com- 
pressible. In the weakly interacting limit, it has been con- 
jectured that the Bose glass goes over into a Lifshits (or 
Anderson) glass [J, Q , where the density profile is that of 
a superposition of the lowest-lying single-particle states, 
which are exponentially localized. Finally, a gas of non- 
interacting bosons will have all the particles occupying 
the lowest of the single-particle states. 

Experimentally, two basic types of disordered poten- 
tial for ultracold atoms have been realized, both by op- 
tical means; speckle and quasiperiodic. The first type of 
disorder uses laser speckle, the two-dimensional diffrac- 
tion pattern of a laser beam passing through a roughened 
plate. Speckle potentials for ultracold atoms were first 
used in Ref . [12]. A quasiperiodic potential in one dimen- 
sion is created simply by superposing two standing waves 
with different wavelength; if the wavelengths differ by an 
irrational factor - the most popular being the golden ra- 



tio - the resulting intensity is a quasiperiodic function of 
the coordinate [l|. In higher dimensions, quasiperiodic 
potentials can be created as the interference pattern of 
a number of beams meeting at judiciously chosen angles 

The speckle potential is intended to achieve the effect 
of a one-dimensional random pattern. This situation has 
been studied extensively in the literature. For a weakly 
interacting (repulsive) system, Albanese and Frolich |14| 
showed that the Hamiltonian admits exponentially local- 
ized solution in the presence of arbitrarily small disorder. 
On the other hand, an interacting ID Bose gas acts as 
a Luttinger Liquid. Using this approach, Giamarchi and 
Schultz |15j showed that for a strong enough interaction 
the system appears again in a localized phase. In order 
to match with the known results in the weakly interact- 
ing limit, they speculated the existence of two localized 
phases. Later on, the existence of those two sides of the 
glass transition has been shown numerically in Ref. [l6[ 
and [13. 

In the quasiperiodic case the situation is less clear. 
While it is established that - for no interaction - the 
system is exponentially localized only when the exter- 
nal potential is strong enough [18[, in the presence of a 



repulsive interaction the exact behavior of the localiza- 
tion is still not known, and only recent results |19| show 
that a quasiperiodic potential can inhibit the diffusion of 
a wavepacket. In this kind of potential the phase dia- 
gram has been studied by Ref. [^, [y| . These works are 
done under the assumption that the secondary lattice is 
small compared to the primary one, and the system can 
be mapped into a Bose Hubbard (BH) model. However, 
very little is known when the secondary lattice strength 
is comparable to the primary one. In this case the BH 
model cannot be used, and one has to address the prob- 
lem in the continuum. Moreover, if there exists an under- 
lying periodic potential to which the disorder is added, 
there is also a Mott insulating phase for strong interac- 
tions if the filling is integer. The latter feature will not 
be discussed much in the present paper because it is not 
correctly described by the approximation that we intend 
to use. 

We see from the above discussion that the correlation 



properties is a key to understanding the behavior of the 
Bose gas; the superfiuid properties is another. Both of 
these can be calculated in the Bogoliubov scheme, where 
the Gross-Pitaevskii equation gives the (quasi-) conden- 
sate density and superfiuid velocity and the Bogoliubov 
equations yield the excitation spectrum and corrections 
to the density. The Bose glass transition in a Gaussian 
correlated potential was explored recently for strong in- 
terparticle potentials by Fontanesi et al. in Ref. Q • The 
Bogoliubov scheme is usually derived starting from the 
assumption that the system is Bose-Einstein condensed 
and there exists a condensate wavefunction $(r), defined 
either as the expectation value of the bosonic field oper- 
ator or as the wavefunction for the single-particle mode 
that is occupied by a macroscopic fraction of the parti- 
cles. Clearly, this assumption would appear to preclude 
the description of both quasicondensate and Bose glass, 
but, as we shall see, it turns out that the Bogoliubov 
scheme can in fact describe those states as well. The ex- 
planation is that BEC is convenient, but not a necessary 
requirement, for deriving the Bogoliubov equations; it 
suffices to make the weaker assumption that the quantum 
fluctuations in the density and gradient of the phase are 
small. Several papers have presented different sketches of 
the derivation of Bogoliubov theory for a quasicondensate 
[2(]| - |25| . In this article we put up what we believe to be 
a complete, simple and consistent derivation of the Bo- 
golioubov equations in a Bose gas using a minimum of 
assumptions, in order to be able to analyze the quasicon- 
densate - Bose glass transition in one dimension. Using 
the same formalism we derive an expression for the super- 
fluid density in order to better understand the transition. 
In this work we consider the experimentally relevant 
case of a ID Bose gas in a quasiperiodic potential and 
examine a wide parameter regime, allowing for a dis- 
cussion of the Bose glass transition and the conjectured 



Lifshits glass. It is important to stress that the Bogoli- 
ubov approximation has to be used with some care in 
one-dimensional systems. This problem has been treated 
in detail for a uniform gas by Lieb and Lininger |26| , 
where it is shown that the Bogoliubov perturbation the- 
ory agrees with the exact answer for g/n < h^/2m. 
Moreover, our numerical scheme is consistent only if the 
fluctuation in the phase and in the density are small, 
and this clearly cannot be the case upon increasing in- 
definitely the interaction within the particles. For this 
reason, due to the intrinsic limit of the Bogoliubov ap- 
proach, we expect our results to be reliable only for weak 
interaction and high density. Within these limits, the ap- 
plication of the Bogoliubov recipe has been proven to be 
theoretically sound, and it is known to give the same cor- 
relation function of the Luttinger liquid theory, as shown 
in Ref. [23| . We wish to compare the different predictions 
for the phase transition obtained from analyzing the be- 
havior of the correlation function and of the superfiuid 
density. 

As we shall see, our analysis suggests the existence of 
a glassy phase for small interparticle interaction. The 
plot of this phase appears as a series of peaks with lit- 
tle overlap. According to the language used in Ref f4|, 
this phase is called "fragmented BEC", or Bose glass; 
in the same reference, the Anderson glass appears when 
the overlap between the peaks is negligible. We maintain 
these definitions through our article. 

This paper is organized as follows. In Section |lT] we 
present a derivation of Bogoliubov theory and in particu- 
lar an expression for a correlation function. In Sec. lIIIl an 
expression for the superfiuid density is derived. Sec. IIVI 
presents numerical results for the superfiuid-Bose glass 
phase transition in a ID Bose gas. Finally, in SecFVlwe 
summarize and conclude. 



II. CORRELATION FUNCTION 

Inspired by the works of Ho and Ma [24I , and Xia and Silbey [23] , we use a path integral formalism, considering a 
system of bosons described by the Euclidean action 



S[r,^A^ j dTdvr{v,T) (-;ia, + A-V2-[/(r)+M~||V'(r,T)|2) V(r,r). 
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where tp is the scalar boson field. Assuming that the quantum fiuctuations of density and gradient of the phase are 
small we expand ip as 



V'(r,T) = e^''("^^"V"o(r) + (5n(r,r) 
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The equation of motion for i/jq is found by means of the variational principle SS/Suq = 0, and the result is the 
Gross-Pitaevskii equation for no [28| , 
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VV'^o(r) + U{r) Vno(r) +gno(r)2 =^^no(r). 
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This amounts to treating no{r) as a scalar field representing the (quasi-) condensate density, and Sn, 6 as perturbations. 
In a three-dimensional system, and in two dimensions at zero temperature, the condensate wavefunction $ would in 
the absence of currents be equal to the square root of the classical density, y/riQ- In principle one could describe a 
state with a current by assuming a slowly varying phase ^o that is also treated classically, but we refrain from doing 
so for convenience. 

Ignoring terms of order higher than 6^ and Sri^ the action becomes 



S — Oo + oi + 02 



(4) 



Here, Sq contains only the classical density tiq and is minimized by the equation of motion for no. Si must vanish for 
uq to be a stationary solution, while 52 is 
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where we have implicitly defined the scalar operators H3 and Hi. In order to find the correlators between Sn and 
we need to find the function G that inverts S, 
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Of course, what we have derived so far is mathematically identical to the well-known Bogoliubov theory, as we now 
show. Introducing the transformation T as 



we obtain 



with 
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The diagonalization of C leads to the Bogoliubov equations for the Bogoliubov amplitudes Uj{r) and Vj{r), 
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Moreover, the Green function for the action S — ihdt + C is already known [29|, and equal to 
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Defining Xj as 
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and correspondingly 



X, - T ( ^i , , 



we can apply the transformation T to Eq. ^ and eventually find 
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with Wn being the Matsubara frequencies, i] a positive infinitesimal, and N{Ej) — l/{e^^j — 1) is the Bose-Einstein 
distribution function. For convenience we have defined the notation 



9' = 0{r',T+) 
6n = (5n(r, r) 
6n' — dn{r' ,T ). 

Furthermore let us define 

M^e' ~ 9. 

Using Eq. ([2]) the one-body correlation function becomes 
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This expression can be evaluated using Wick's theorem. To lowest order in 6n and one finds [24 

,dn 

-3 ((AS) 



, (5n 



"0 / V "^0 , 



no 

2 



,dn 



no /tg 

so that, eventually, to second order J23 l. 
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Since both 5n and are small, the expression between square brackets can be thought as a first order expansion of 
an exponential. Using Eqs. (IT5t and P^ the following expression is obtained, 
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where we useed the shorthand Nj — N{Ej). The quantity il{r,r') is purely imaginary, and equal to 
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If the excitation energies and the condensate wavefunction are real, then the Bogoliubov amplitudes Uj and Vj can 
also be chosen as real. The imaginary part in Eq. (I22p therefore vanishes, and the resulting expression coincides with 
that found by Mora and Castin using a different formalism [2J|. In particular, for the uniform case the expressions 
for the Bogoliubov functions are known, leading to 



In one dimension, Eq. ()2ip is an expression for the correlation in a Bose system which is ultraviolet and infrared 
convergent. Even at T = the sum has a finite value without assuming a cutoff or a modification in the interparticle 
potential. 

III. SUPERFLUID DENSITY 



According to the two- fluid model j30l |31| , the mass density of a quantum fluid can be divided into a superfluid part 
Psiy) and a normal one Pn(r), the total mass current being J(r) = p„(r)ws(r) + p„(r)w„(r). While the superfluid 
density is in general different from the condensate density, the superfluid velocity is the condensate velocity. In 
particular, upon imposing a phase twist on the condensate wavefunction, the superfluid part will be proportional to 
the additional kinetic energy. 

It is in this sense that the superfluid density can be defined as a response to a twist of the order parameter [32, ISJ] , 
by means of rewriting Eq. ([2|) as 

^(r) = e^^WV7^fe^''o- + i^-l^), (24) 

V 2 no 8 ng / 

with ko = eo/i, h being the length of the system in the direction of the unit vector eg, and G a small twist angle. 
For convenience, let us take the order parameter normalized to unity; the superfluid density - in the direction of cq - 
is then defined by the thermodynamic limit of 

2L2to2jy 
P^ = f,2Q2v [F''if^,T)-F\^,T)] (25) 

The substitution (p4)) results in the twisted action 

f h^ 

S^ = S+ drdT \k?,no(r,T)-2ikaSn(r,T)\/i9(r,T)] 

J 2m '- ■' 

= S + ko fdrViT), (26) 



where in the second hne we ignored higher order terms in the fluctuating fields Sn and 0. If the system is large enough 
{L 3> 1), then V{t) can be seen as a perturbation; moreover, let us assume that the symmetry of the problem imposes 
that the odd terms in fco vanish in F®, since the cases with ko and —ko lead to the same physical situations. The 
new free energy can be computed using the linked cluster theorem (see, for example [3J]) 
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Here, the notation (. . .)connccted stands, as usual, for a diagrammatic expansion where only connected diagrams are 
retained. The superfluid density can be found in terms of the Green functions (J15l) , defining 
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so that, by partial integration of Eq. ([27| . 
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The evaluation of the averaged operators is done by summing over the Matsubara frequencies in expressions like 
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with Xi and Xi defined in Eq. P^ and (|14l) . In particular, for the uniform case, all the terms proportional to Vy^no 
vanish; moreover, since the expressions for the Bogoliubov functions are known analytically, it is possible to see that 



only the last term of ((30)) survives, each term of the sum between square parenthesis giving a contribution of fc^/2. 
In one dimension one obtains 
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which is the well known Landau result 13 ill 
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IV. BOSE GLASS TRANSITION 



(31) 



We now apply our expressions for the correlation func- 
tion and the superfluid density to a system at T = 0. The 
external potential is quasiperiodic, obtained as a sum of 
two potentials whose periods are incommensurate with 
each other, 



'27r 

U(x) = V\ cos ( — - X 
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(32) 



For our specific realization we have chosen A to be the 
golden ratio, approximated as a fraction of two consec- 
utive Fibonacci numbers. In the following we shall con- 
sider two systems with length L — 377 d and L = 89d; 
the value of A is, respectively, A = 377/233 and A = 
89/55. 

For V2 = 0, upon increasing the value of Vi, the system 
ceases to be superfluid and enters a Mott insulating phase 
[lOl I35I l3q |. As stated in the introduction, this phase 
cannot be seen using the Bogoliubov ansatz we employ 
for the excitations, and our method should not be applied 
when the physics of the system is affected by the presence 
of the Mott phase. However, when the second lattice is 
turned on a new quantum phase, the Bose glass, becomes 
possible. According to Ref. [l^, when V2 is greater than 
the interaction energy, the Mott phase disappears from 
the phase diagram of the gas, and the only insulating 
phase is the glass phase. We believe that in this situa- 
tion of strong disorder our approach can be quantitatively 
correct, and for definiteness we have chosen the weights 
of the two potentials to be equal (Vi — V2 = V). 

Taking d as the unit of length, the Gross-Pitaevskii 
equation can be rewritten as 



l^+gNmx)\' + U{x) 



^x)^^inx), (33) 



where the energy is measured in units of Ed = h^ /md^, 
and the norm of the quasicondensate wavefunction $ is 
set to 1. For each value of the parameters V and 5, 
we have found the ground state of Eq. ()33|) using an 
imaginary time evolution with the split-operator method, 
along with periodic boundary conditions. Using the 
ground-state wavefunction the Bogoliubov excitations 
arc found by means of a direct diagonalization of the 
Bogoliubov equations (jlip . In all our calculations the 
Laplacian term is represented using the Fourier trans- 
form, while the derivative in the expression for superfluid 
density (|29p is approximated using a finite difference ex- 
pression. 




FIG. 1: (Color online) Density profiles for some values of gn 
for V = 5 Ed (top panel) and V — 9 Ed (bottom panel); <; n is 
expressed in units of Ed ~ T:? /mS and x in units of d. Upon 
decreasing the interaction strength the condensate breaks up 
in sets of spikes with little overlap. 



The Bogoliubov analisis is expected to be relevant for 
values of the interparticle interaction g/n <^ 1. This 
means that our work is realistic for 



g n 't^ n 



(34) 



i.e., the approximation we are using starts to be mean- 
ingful for systems that contain a few particles per site. 

We first solve for the ground state in a system with 
L — 377 d, using a numerical grid of 3016 points. Choos- 
ing a mean density n — N/L — 4, we show in Fig. ([T]) 
no{x) = 7V|$(a;)p for the two cases V — 5 Ed and 
V = 9 Ed- Upon decreasing g, the condensate density 
develops dips that become more and more pronounced. 
Indeed, in some regions the condensate seems to be bro- 
ken up into several pieces that hardly overlap. As we 
shall see, this has consequences for the behavior of the 
correlation function. 

Using the same mean density n = N/L = 4, Fig. ([2]) 
plots the exponent of the correlation function 

lngi{0,x) = ln(V'(0)i/^(a;)) -1/2 In no{0)no{x) . (35) 
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FIG. 2: (Color online) Exponent of the correlation function 
((35} for the same values oi gn as in Fig [l] for V = ^Ed (top 
panel) and V = %Ed (bottom panel); gn is expressed in units 
of Ed = Y? jmS and x in units of d. The insets shows the 
same plot with a log scale on the horizontal axis. In the main 
picture it is seen that for lower gn the values of In g\ psp 
follow a linear behavior. For higher gn the plot appears to 
depend logarithmically on x. 



While for higher g n the values of In gi follow a logarith- 
mic behavior (insets of Fig. ^ ) , leading to a power law 
decay in the correlation function, for lower g n the func- 
tion shows a linear fall, therefore giving an exponential 
decay. Since our system is finite the transition to a dif- 
ferent decay behavior is gradual, but as a quantitative 
measure we record the point where a linear fit for In gi 
gives a smaller sum of the residuals than a logarithmic 
fit. The full line in Fig. ^ shows the critical value of V 
as a function of the interaction strength, according to an 
interpolation of In gi between x = 10 d and x — 150 d. 

The change in the decay behavior becomes evident at 
the values of g when the quasicondensate seems to be 
broken up into a set of spikes with little overlap. In this 
situation the lowest energy excited modes have a phase 
flip character, since it costs little energy to change the 
relative phase of two barely overlapping zones; the en- 
ergy of the lowest excitation decreases upon decreasing 
g. Eventually, this value becomes so low that reaches 
the limit of our numerical precision (^ 10^^^ Ed). At 
this point, we may say that the condensate has become 
completely disjoint. This boundary is indicated by the 
dashed line in Fig. ([3]). It is seen that it occurs close 
to the Bose glass transition, but it does not coincide 
with the transition found using the other criteria. In the 
lower left corner of the phase diagram, for small values 
of the interparticle potential, detecting the transition can 
be numerically challenging, since at the point where the 
correlation function changes its behavior the excitation 
spectrum is already at the boundaries of the numerical 



FIG. 3: (Color online) Phase diagram of the Bose glass tran- 
sition. The open circles correspond to a state with non- 
vanishing superfluid component, while for the red dots the 
system has only a normal part. The filled line indicates when 
the linear interpolation of In gi (Eq. (|35|) ) gives a smaller er- 
ror than a logarithmic fit, i.e., at the left side of the filled line 
the correlation function decays exponentially. The dashed line 
indicates the point at which the smallest excitation energies 
cannot be resolved numerically. The inset shows the sum of 
the residuals of the fits (in arbitrary units) as a function of 
gn, for V — 9 Ed'- the triangles refer to a linear regression, 
the diamonds to a logarithmic interpolation. V and gn are 
expressed in units of Ed ~ h^/mcP. 



precision. 

Figure ^ also shows when the superfluid density, as 
given by Eq. (1^^ . drops to zero (red points). Apparently, 
the superfluid formula systematically overestimates the 
critical interaction strength for the transition to an ex- 
ponential decay. Indeed, Eq. (P^)) gives the superfluid 
density only in the thermodynamic limit and it might be 
that the finite size of the system is part of the expla- 
nation for the disagreement. However, we believe there 
is another reason for this discrepancy: The way we ap- 
plied the phase twist to the condensate part in Eq. ([21]) . 
For a non-uniform system, a phase twist in the boundary 
conditions should give a velocity that depends on the co- 
ordinate, as a function of the condensate density in each 
point. The corrections to (P^ should be more relevant 
for small g, when the variations in density are bigger. 

The interpretation of the phase diagram in Fig. ^ 
is the following: at the left side of the line there is an 
exponentially decaying correlation function, along with 
a vanishing superfiuidity; we believe that at this point 
our Bogoliubov scheme detects the Bose glass phase, as 
described previously in Sec. HI The Bose glass phase is 
seen to disappear completely below a finite value for the 
disorder strength V (approximately V^ = 3), in contrast 
to the findings of Ref. 7]. This is the main difference 
with the results of Fontanesi et al.. Besides finite size 
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FIG. 4; (Color online) Surface plot of the superfiuid density 
(|29|l (normalized to the total density n) as a function of V 
and g n, for L = S9d. The dashed line is the boundary that 
we estimate for the validity of our numerical calculations (see 
text). The inset shows the inverse participation ratio P p6|) 
of the quasicondensate wavefunction. Both V and gn are 
measured in units of Ed = Ti' jmS . 



effects, this behavior seems to be characteristic of the 
quasiperiodic potential. As stated in Sec. U for vanishing 
interaction the system is localized only when the external 
potential strength is higher than a critical value. When a 
repulsive interaction is considered, this property survives 
in the glassy phase. 

In Fig. m we show the behavior of ps in an enlarged 
phase diagram. Because of the computational burden, 
this phase diagram was calculated for a smaller system 
with L = 89(i and a grid of 1424 points, and therefore 
the Bose glass boundary is slightly shifted compared to 
that in Fig. [31 Monte Carlo simulations [la, Hzt have 
previously shown a reentrant phase transition upon in- 
creasing g. While we do not see this re-entrance, the 
behavior of the superfluidity is suggestive, as it shows a 
non-monotonous dependence of the normal phase with 
respect to the interaction and the external potential 
strength. Indeed, for a quasiperiodic potential the pa- 
per Ref. fq] shows that the re-entrant phase is limited to 
a small portion of the phase space, and it is compatible 
with a superfiuid phase that extends to higher interaction 
strengths. Within the limit of the Bogoliubov approxi- 
mation, our work agrees with this conclusion. 

Finally, Fig ([5]) and the inset of Fig. ([4]) shows the 
inverse participation ratio of the quasicondensate wave- 
function $, 



L /dx|$(a;)|4 



(36) 



for both of the cases L = 377 d and L = 89 d. Since 
the computation of P is much faster than the diagonal- 



FIG. 5: (Color online) Inverse participation ratio P (|36p for 
L — 377 d as a function of g n, for some values of the external 
potential strength V. The inset shows a surface plot of P as a 
function of V and g n. All the energies are measured in units 
of Sd ^h^/m<f. 



ization of the Bogoliubov equations, Fig. ^ presents 
the behavior of the system for a finer resolution in the 
values oi gn and V than Fig. ([3]). As we can sec, the 
plots do not show any distinct transition. The glassy 
phase does not appear to be related to a radical change 
in the density profile. Moreover, as stated in Sec. HI a 
new phase, the Lifshits glass, has been conjectured to 
exist for weak interparticle potential, where the peaks in 
the density profile do not overlap. We argue that such 
a phase should have a distinctive signal in the quantity 
P. However, the plots of the inverse participation ratio 
do not show only a smooth decrease. We conclude that 
- within our approximation - the Lifshits glass does not 
seem to exist as a phase in its own right, but only as the 
limit of vanishing interaction. 



V. CONCLUSIONS 

Using a path integral formalism we have obtained - 
within the Bogoliubov approximation - expressions for 
the correlation function and the superfiuid density in a 
Bose gas for an arbitrary external potential, at zero or 
finite temperature. Applying these expressions to a one- 
dimensional system at T = 0, a quasi-periodic external 
potential was seen to cause a transition to a phase where 
the correlation decays exponentially and the superfiuid 
density vanishes. We believe this is the Bose glass phase. 
A comparison with the experiment of Fallani et at. [1] is 
not straightforward, because they employed a quasiperi- 
odic potential with different relative weights [i.e., Vi ^ V2 
in Eq. ([5^ ]. and the proximity to the Mott insulating 
phase in the experiment makes the Bogoliubov approx- 
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imation questionable. However, they spotted signs of a 
Bose glass transition when the peak height of the exter- 
nal potential varies between 2V fv 16 Ej, and 2V ~ 18 Ed, 
with gn fa 1 Ed, which is where the transition takes place 
in our study for the bigger system (L — 377 d). 

We have not found a perfect match between the results 
for the superfluid density and the correlation behavior for 
the location of the Bose glass transition. We would like 
to point out that the formula for the superfluid density 
([^ assumes that the velocity Vg of the condensate does 
not depend on the position, and in this sense it is an 
approximation even within the Bogoliubov approach. For 
this reason we plan to improve the expression for ps in 
the future. 

Finally, we notice that our expression for the correla- 



tion function and the superfluid density are valid also at 
a finite temperature, and they can be used to study the 
Bose glass transition for T 7^ 0. This relevant issue will 
be considered in a future work. 
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